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In this paper we discuss the well known Kardar Parisi Zhang (KPZ) equation driven by tem- 
porally correlated noise. We use a self consistent approach to derive the scaling exponents of this 
system. We also draw general conclusions about the behavior of the dynamic structure factor 
& q (t) as a function of time. The approach we use here generalizes the well known self consistent 
expansion (SCE) that was used successfully in the case of the KPZ equation driven by white noise, 
but unlike SCE, it is not based on a Fokker-Planck form of the KPZ equation, but rather on its 
Langevin form. A comparison to two other analytical methods, as well as to the only numerical 
study of this problem is made, and a need for an updated extensive numerical study is identified. 
We also show that a generalization of this method to any spatio-temporal correlations in the noise 
is possible, and two examples of this kind are considered. 
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I. INTRODUCTION 



Nonequilibrium surface growth processes often exhibit a phenomenon called kinetic rough- 
ening, where the surface develops a self-affine morphology jlj . Much attention has been given 
to a special class of models (ballistic deposition, Eden, or polynucleation growth), which are 
described by the Kardar-Parisi-Zhang (KPZ) equation Q] 

dh(f,t) A,„,, 2 ,_ . 

—±±l = uw 2 h + - iyhf + n (f, t) , (i) 

where h (r, t) is the local height of the surface above a <i-dimensional substrate in a {d + 1) - 
dimensional space, A characterizes the tilt dependence of the growth velocity, v is an effective 
surface tension, and n (r, t) is a noise term. 

Solutions of eq. (JTJ) exhibit scaling behavior. The simplest quantity to investigate is the 
surface width W (L,t) that scales as (see ref. 0) 

W(L,t) = ±(£[h(r,t)-h(t)] 2 y =L a g(J^), (2) 

where h (t) is the mean height of the interface at time t, a is the roughness exponent of the 
interface and z is the dynamic exponent that describes the scaling of the relaxation time 
with L - which is the size of the system. The brackets (• • •) denotes noise averaging. The 
scaling function g (u) behaves like g (u) ~ u 13 (where f3 is the growth exponent) for small 
u's (i.e. for t L z ) and like a constant (i.e. g (it) ~ const) for large it's (i.e. for t ^> L z ). 
It is easily verified from eq. (0) that (3 = a/z. The scaling exponents a and z describe the 
asymptotic behavior of the growing interface in the hydrodynamic limit. 

The KPZ equation with uncorrelated noise has been well studied. For the one dimensional 
case one can easily obtain ex^ct results of a = 1/2 and z = 3/2 by mapping the KPZ equation 

by using the Fokker Planck equation associated with the 
However, for higher dimensions (d > 1) there are no 
exact results and the critical exponents have been evaluated numerically or obtained using 
various analytical methods (for a review see 



into the Burgers equation (2J or dy 
Langevin form given by eq. (JIJ l| 



The noise in the KPZ equation is a result of a physical process. As such it must be 
correlated in space and in time. If the correlations in space and time are short ranged it may 
be expected that the long distance and the long time behavior of the system characterized 
by the exponents a and z are those obtained in the case of uncorrelated noise. There may 
be, however, situations in which the decay of correlations in the noise is algebraic. 



Indeed, in some experimental situations the measured scaling exponents are larger than 
the values predicted by KPZ 0,0. A possible explanation of such a departure from KPZ 
behavior may be long-range correlations in the noise. Such experimental results serve as 
serve as a motivation for the study of systems with correlated noise in spite of the fact that 
direct evidence for long range correlations in the noise is usually lacking. 

Many studies of growth models with noise that is algebraically correlated in space but 
uncorrelated in time described by 

(V(r,t)) = 0, (3) 

and 

(r]{r,t)r}{f>,t)) = 2D \f - 7\ 2p ~ d 6 (t - t') , (4) 
have been published in the last decade. These include discrete one-dimensional models 



:bd 0,00, sos 00, 



], and direct (discrete) integration of the KPZ equation |5| ). Many 
researchers studied the KPZ equation with such noise and 



obtained different predictions. In spite of the differences in the predicted values of the critical 
exponents, a common picture seems to result from all methods, namely: for small p's the 
critical exponents are the same as for the case of uncorrelated noise. Then, for p's above a 
certain critical value p c the exponents become p-dependent. 

In sharp contrast to the variety of numerical results and theoretical predictions for the 
critical exponents of the KPZ equation with spatially correlated noise, only few results are 
available for the KPZ equation with temporally correlated noise - not to mention noise that 
is both spatially and temporally correlated. Similar to eqs. ©-(jlj), temporally correlated 
noise with zero mean can be described by 

( V (r, t) V (r"*, t)) = 2D (f - ?) \t - t'^ 1 , (5) 

where <ft characterizes the decay of the correlations over time (it is assumed that < 1/2 or 
otherwise the correlations does not decay, but rather increases with time). 

The first theoretical prediction of the critical exponents of KPZ in the presence of this type 
of noise is due to Medina et al. [13] that used Dynamic Renormalization Group (DRG) anal- 
ysis to study this problem. They solved the DRG equations numerically in one- dimension, 
for the case where Dq is a short range function, and found out, just like for spatially cor- 
related noise, that for small enough 0's the correlations are irrelevant. They claim that 
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for > 0.167 the correlations become relevant, and the roughness exponent can be fitted 
numerically to 

o^DRG (0) — 1-690 + 0.22. (6) 

The dynamic exponent can then be obtained using the scaling relation 

2a DRG ((f)) + 1 

Z DRG (0) = 1+ ■ (7) 



These predictions have been checked numerically by Lam et al. [17( using the Ballistic 
Deposition model. They found sensible agreement between the DRG prediction and the nu- 
merical values they obtained. However, substantial deviations were found, centered around 
the expected threshold point 0o = 0.167. Thus, the authors believe that these discrepancies 
are due to a crossover effect in the simulation and not due to any approximation in the DRG 
calculation. 

Apart from the above mentioned DRG result, there is only one more result for KPZ with 
temporally correlated noise due to Ma and Ma |Tj| who used a Flory-like Scaling Approach 
(SA), originally suggested in the white-noise KPZ context, by Hentschel and Family jl^ . 
Ma and Ma obtained the following strong-coupling roughness exponent 

/ n 2 + 40 
^ (0) = 20 + rf + 3' (8) 

and the following dynamic exponent 

^ (0) = 20T7+T (9) 

These values are said to describe the strong-coupling scaling exponents for all values of the 
parameter 0, and for every dimension d. Actually, it is easily verified that these expression 
reduce to the well-know white-noise KPZ results in one-dimension when = 0. 

This prediction for the critical exponents is obviously different from the previous one- 
dimensional DRG result in two respects. First, Ma and Ma do not predict that for small 
enough 0's the temporal correlations are irrelevant, so obviously they rule out the threshold 
value of 0, 0o- Second, for > 0.167 the two approaches yield different numerical values for 
the scaling exponents. 

This situation, where only two theoretical predictions are available for the KPZ problem 
in the presence of temporally correlated noise, especially when one of them (DRG) is a one- 
dimensional result, certainly calls for a clarification of this issue. This problem is further 



complicated by the fact that only one numerical study [171 ]. and only in one dimension, is 
available. 

At this point it is interesting to mention another result for the KPZ equation in the pres- 
ence of noise with special mixed spatio-temporal correlations (non seperable noise correlator 
D(q, uj)). This is a case where in contrast to systems where the noise is only suspected to be 
of long range, here long range correlations in the noise follow from direct physical arguments. 
This problem has been studied both numerically and analytically by Li et al. ^| with good 
agreement between the analytical and numerical values. Since we deal with this problem in 
section VI we will not discuss it further now. 

In this paper we develop a self-consistent approach to deal with nonlinear Langevin 
equations, such as KPZ, with temporally correlated noise. Actually, as will be seen in 
section VI, this approach can be easily generalized to spatio-temporally correlated noise. 
We begin with a brief derivation of the scaling exponents of the linear theory (also known 
as the Edwards- Wilkinson equation) in the presence of temporally correlated noise. Then, 
the full time-dependent two-point function for the linear problem is derived. This result 
will serve as a reference for the more general nonlinear discussion. In section III concepts 
emanating from a previous self-consistent Fokker-Planck expansion to the KPZ equation 
are reviewed. In section IV the time-dependant self-consistent approach is established. It 
is shown that analysis of the time-dependant self-consistent equation in the limit of short 
times and long times yields two static equations that are an interesting generalization of the 
former self-consistent Fokker-Planck expansion. 

In section V a detailed asymptotic solution of the self-consistent equations is obtained. 
In this section, we derive the different possible phases and their corresponding scaling ex- 
ponents. Special attention is given to the results in one dimension. Section VI generalizes 
the previous results to the case of noise with arbitrary spatio-temporal correlations, and two 
elaborated examples are given. At the end, in section VII a brief summary of the results 
obtained in this paper is presented. 

II. THE LINEAR THEORY - THE EDWARDS- WILKINSON EQUATION 

At the beginning of this paper we would like to discuss first the Linear theory (i.e. the 
KPZ equation with its coupling constant set to zero - A = 0), namely the Edwards- Wilkinson 
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(EW) equation [20], with temporally correlated noise. The Edwards- Wilkinson equation is 

dh 

-^(r,t) = vV 2 h + V (r,t), (10) 

As mentioned above, in this paper we discuss temporally correlated noise characterized by 

( V (f,t)r](f^,t')) = 2D 5(f-f')\t-tf' t " 1 , (11) 

where the case of uncorrelated noise corresponds to the limit = 0. 

The interface that grows under these conditions is known to be self-affine, which means 
that if the spatial coordinates are scaled by a factor of b (i.e. r — > r 1 = br) then if we 
perform the transformations t —>■ if = b z t and h —>■ hf = b a h (with the appropriate scaling 
exponents a-the roughness exponent, and z-the dynamic exponent) as well, the statistical 
properties of the surface are left invariant. Since the growth equation is linear, following 
ref. it is possible to extract the scaling exponents by scaling r, t and h in the equation 



according to the above-mentioned transformation. But first, we have to realize 
this transformation the noise term scales like r\ — ► rf — foW 2 '^ 1 ) -0 ')/ 2 ^ (see ref. 
this we can plug it back into the EW equation and we get 



that under 
) . Using 



i] 



f)h' 

b a - z jjj {f, t') = b a - 2 vV 2 ti + b^-^-^r]' t') (12) 

Now, imposing the requirement that eq. (|10|) remains invariant under this scaling transfor- 
mation, namely requiring that both equations (eqs. (fTUj) and should be exactly the 
same, we get 

z = 2 and a = (40 + 2 - d)/2. (13) 

(this gives the roughness exponent as long as the resulting a is positive, otherwise the surface 
is flat). It is easily seen that this result reduces to the standard EW exponents (i.e. for the 
EW equation with uncorrelated noise) in the limit of = 0. 

This simple result shows that temporally correlated noise tends to make the surface 
rougher (a bigger roughness exponent a implies a rougher surface). 

The information extracted so far regrading the EW equation in the presence of tem- 
porally correlated noise could have been satisfactory. However, because we are interested 
in obtaining the exponents of the nonlinear theory as well, we would like to gain as much 
insight into the behavior of the linear problem, so that it might help us when dealing with 
the KPZ nonlinearity. For example, because of the linear character of eq. (jlOj) . we can 



obtain the scaling form (and recover the exponents) by solving the growth equation exactly. 
Fourier transforming eq. (fTU|) in space and time we obtain 



h 



(14) 



iuj + vq 2 

where rj qLJ is the Fourier transform of rj (f, t). Thus, using the Fourier transform of eq. ([lip, 
we obtain the dynamical structure factor (or the two point correlation function) 

-2c/> 



qui 



2 n 4 



(15) 



By Fourier transforming back we get 

(Jt) = {h q (0)h- q {t)) 



D 



u l+2<f> cog (ft. 



7 (m 2 t 



(16) 



Here f(u) is a scaling function that can be written explicitly as 



f(u) 



COS (7T<; 
7T 



; -e^dy = cosh (w) - 



Z/ 2 + l 



T (2 + 20) 



1-^2 



V + 1,0 + 



^ ■ (17) 



where T(x) is just Euler's Gamma function, and 1F2 is a generalized hypergeometric function. 
The function f(u) is also plotted in Fig. ^ As can be seen in the figure, the scaling function 




FIG. 1: The scaling function f(u) (0 = 1/4 was taken for this illustration). One can see an 
exponential-like decay for small u's, and a power law decay for large u's. 

behaves like a constant for small u's (this corresponds to short times, that is for vq 2 t <C 1). 
At the other extreme, i.e. for large u's, this function decays algebraically. In order to be 
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sure of this power law tail, and to obtain its exact shape we calculated the leading behaviors 
for small and large it's and obtained 



-u 




1+2 * H u < 1 



/(«)~< „*-i L *<W)(§-*) \ , (18) 

J? 1 I M > 1 

so that /(w) ~ ii~( 1_2< w for large it's. 

Naturally, the scaling exponents can be recovered easily from $ 9 (i). Since $ g (i) depends 
on time only through the combination uq 2 t we identify the dynamic exponent as the power 
of q in this scaling form, so that here z — 2. In addition, it can be seen that for small g's, 
&q{t) ~ q~ 2 ~ 4 ^- Thus we identify the exponent T = 2 + 40 that can be translated into the 
roughness exponent via the relation a = (T — d)/2 (see eq. ()46|) below), so that we recover 
a = (40 + 2 - d)/2. 

The results obtained in this section will serve us later. First, it might be interesting to 
compare these results with the results obtained for the nonlinear theory (for example, in 
the weak-coupling regime of the KPZ equation). Second, we will use the scaling function 
of the linear theory as an ansatz for the integral equation that will determine the scaling 
exponents of the strong-coupling phase of the nonlinear theory. 



III. NONLINEAR THEORY - THE KPZ EQUATION 

We proceed now to the much harder nonlinear case that poses many technical difficulties 
already in the uncorrelated case. 

The method we present in the following section is based on the same general ideas as 
;he self consistent expansion (SCE) used for systems with noise that is uncorrelated in time 

inn 

9|, |2a, |26(. Namely, an expansion around an optimal linear system. The SCE is based on 
constructing a Fokker-Planck equation for the probability distribution of the height function. 
This step is based on the fact that the noise is not correlated in time. The self consistent 
expansion is formulated in terms of the steady state structure factor (or two-point function), 
4>q — (hqh-q) s an< ^ ^ s corresponding steady state decay rate that describes the rate of decay 
of a disturbance of wave vector q in steady state 

-i SS°(h q (t)h^(o))dt 

^ - • (19) 
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The linear model around which the expansion is constructed is chosen to yield the (un- 
known) <p q and uj q that appear in it as parameters. An evaluation of <p q and u q as an 
expansion around that linear model leads to the coupled equations 

<Pq = 4>q + C q {(j) p , Up} , (20) 

and 

U q = UJ q + d q {<f) p , U p } . (21) 

Within this framework, the structure factor and decay rate are obtained by solving the 
coupled non-linear integral equations c q {<p p , u p } = and d q {4> p , u p } = 0. In contrast to 
other expansions, the full correction, in a given order of the expansion, of the relevant 
physical quantities, is really small. In fact, it is chosen to be zero. 



IV. DERIVATION OF THE TIME DEPENDENT SELF-CONSISTENT AP- 
PROACH 

In this work (following ref. [22']) we obtain the dynamical structure factor = 
(h q (0) h_ q (t)) , using the same idea of a self consistent expansion. Here too, the average 
(■ • -) s denotes steady-state averaging, where h q (0) is measured in steady state at time t = 
and then h- q (t) is measured at some later time t (also in steady state). The dynamical 
structure factor $ g (t) normalized by 4> q = $ 9 (0) (i.e. the static structure factor) is thus a 
measure of the persistence in steady state of disturbances with wave vector q. Because the 
noise is correlated in time, we cannot use the Fokker-Planck approach, but as seen in refs. 



22 



23| such an approach lends itself as an alternative to the Fokker-Planck approach even 
when it is available. 

Our starting point is the field equation for h qLU (the Fourier transform in time and space 
of h (r, £)) obtained by Fourier transforming eq. (JTj) 

iuihqu + iSqh qu > + C q e m he a h mT = r] quJ , (22) 

where u q = uq 2 , C q £ m = -jf-j^^ • w^g,i+m ( L,<T+T, T being an assumed periodicity in time to be 
taken eventually to infinity, Q is the volume of the system (to be taken to infinity as well) and 
the noise correlations are (^^-g-w) = 2-D (q) u~ 2 ^. (Note, that in the d + 1 dimensional 



9 



space (including time), the noise is quenched disorder!). In the Chapman- Enskog spirit (as 
done in refs. |2ll. |22|) the equation is written in the form 



(iu + u q ) h qa} - r]° quJ 



+ A 



+ A 2 [{v q - uu q ) h qul ] = 0, (23) 



where A is going to be taken as 1 but is used at present as an indicator to show the construc- 
tion of the perturbation expansion as an expansion in A. The noise is split into two terms 
Vqui — Vqw + r 7gu; sucn that (j]quV-q-u/ = Dqu and the correct (i.e. the Fourier transform 
in time of the "dynamical structure factor" $ g (t)) is given by § quJ = w2 j^2 ■ This choice 
implies that ignoring the A and A 2 terms in eq. ([23)1. we still obtain from a linear equation 
the correct In contrast to the case of short range correlated noise where u q is defined 
by eq. ()19|) . we must employ here a more general definition. The reason is that the power 
law found to describe the tail of & q (t) for long times renders the expression on the right 
hand side of eq. (fTT)|) infinite. Therefore, our definition of uo q is based on the assumption of 
a scaling form of - namely 

* q (t) = </> q f(u q t), (24) 

(it can be easily verified that the dynamical structure factor of the linear theory given by 
eq. (fTBj) indeed obeys this scaling law). The "decay rate", u q , is defined as that parameter 
that will make eq. ()24|) a good approximation for small g's and over the whole time range. 

Now, Eq. (J2l!j) enables to obtain h qu} explicitly to second order in A. The expression for 
h qw is multiplied into its complex conjugate and only terms up to second order in A are 
retained. At the end the expressions are averaged over r) qiJ and we get 



(u 2 + u 2 ) <S> quj = D quJ + 2A 2 C, fc C H -A^ 

£,m,a,r 



+ A 2 [2D° (q) u~ 24> -D q J- 2A 2 (v q - u q ) u q % 



+ 4A 2 ^2 C q e m Ce q - m & qoj & r 



£,rn,a,T 



■ mr 

lO + Uj£ 



+ 4A 2 V" C q ^i^ m Ci^ qm Q mT Q quj — ; — ■ — — . (25) 

o —la + Up 

Now A is set to be 1. The result is an equation of the form Q quJ = <& quJ + e {&e a }- Equating 
e {$£0-} t° zero yields 



U) 2 + u\ + 2 (V q - UJ q ) U q 



$ ?a , - 2 J2 \ C lirn\ $fo$mr + 
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+4 ^2 C q e m Ct qm $ qu! $r 



—iuj + uj a 



+ 



iuj + uj a 



2D (q)uj- 2 ^. 



lO + UJ, —lO + Ld£ 

We divide the last equation by (uj 2 + uj^j and using the definition C q e m = Al 'i- iS ^ l+m 
well as letting fl and T tend to infinity) we obtain 

<!''( do A 12 



(26) 



as 



1 + 2 



v Z/g - UJg) UJg 

UJ 2 + oj 2 



+ 4 



/ 



_d^da 
(2^27' 

-20 



e,q-eslq,q- 



1 1 (2n) d 2n 

Qqu^q-^u-a 



UJ 2 + tf2 



+ 



[ioj + u; g ] [ia + t^] [— icr + toe] [—iuj + a> g ] 



uj 2 + uj 2 



2D° (q) . 



(27) 



The last equation is the basic equation for our following discussion. We consider first the 
small uj behavior (more specifically uj/uj q <C 1) that corresponds to the long time decay of 
the time dependent structure factor. 



A. Long time decay of the structure factor 



The first small uj simplification is obtained by neglecting uj/uj q . This yields 



1 + 2 



'Vq ~ V q ) UJ q 



OJi 



n 1 f d d i da . . |2 ^ ^ 

$ ow - 2— / -r \At a -t\ ®( a $ -f 

uj 2 J (27r) d 27T 1 l ' q " q 



®qu, f d d £ da uj, 



2^° (g) , -20 
t: UJ 



uJi 



Fourier transforming back from frequency domain uj to real time t, we obtain 

2 



1 + 2 



Vq - Uq) ^q 
UJ 2 



•.<«> - ^/ 



+ 



q- (2* 
8 r d d i 



(28) 



UJn 



(2n 
D° (q) 



l\Ae, q -e\ *t (*) %-i (!) 

/■no 

dt'e-^'^e (t')$ q (t-f) 

o 

r (i-20) ; (29) 



1 A l , )q _ £ Ag / _g 



T (20) COS (71(f)) UJ 2 

where on the right hand side we have written only the leading large-t behavior. 

This result suggests that in the long-time limit, the time-dependent two-point function 
has an algebraic decay of the general form 



% (t) ~ A^<Pg (Ujgty , 



(30) 



11 



where A°° is a numerical constant, <p q is the steady-state two-point function, and 7 is an 
exponent that will be determined later. 

Equipped with the last result we can see that the first integral on the left hand side of 
eq. (|2fjjl is negligible compared to the other terms on that side in the long-time limit. The 
reason is that this integral decays as t~ 27 , while the other terms decay as t~ 7 , making that 
integral subdominant for large t's. 

Next, using this simplification as well as the scaling form (I24[) . we analyze eq. (|2~9j) for 
small q's (i.e. in the large scale limit) in the spirit of refs. 0, b^ . In order to achieve that, 
we break up the integral into the sum of two contributions corresponding to domains of i 
integration, with high and low momentum. When performing this under the assumption of 
long-times (i.e. u q t ^> 1) we obtain the following equation 



2u q - UJ q 

UJ„ 



<t>q 



(2tt) ( 



f(uj q t) 



'/.) 



D iq z + / <flA. 



',q-eA q 



({/} 



D° (q) (ujq 



-l-2<t> 



(31) 



A°°T (2(f)) cos (ircj)) 

where /, inside the curly brackets, is just the scaling function, q is the upper cut-off of the 



small 



region 



region, and Di is a constant that comes from the contribution of the large 
of the first integral (see ref. j3] section VI, where such an estimation of the contribution 
of large momenta is also employed). In addition, we used the following notation for the 
integral F\ 

Fi ({/} 



00 _j£e 

e 





-(1-20) 



i~ e f (x) dx. 



(32) 



We conclude that (t) ~ A°°0 q (uj q t) (L ' m (i.e. 7 = 1- 20). And eq. flUTJ) can be 
re-written as a time-independent equation relating the 0's and the a>'s 



2u q + D ig 2 -UJ q + 
UJq [ {2tT) 

D» (q) (UJ^ 



d d iA f „_,A 



e,q-esiq,e-q 



<Pq-l 
Uq-t 



Fi{{f}, 



(33) 



T (20) COS (7T0) ' 

It is interesting to compare the result above for <& q (t) with the decay in the case where 



the noise is not correlated in time (namely, when 
behavior of is given by 



0) 



22j. In that case the long time 



$ 9 (t) OC (Ujqt) 2 * exp - (Ujqt) 



(34) 
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(i.e. a stretched exponential). 

The limit as <fi tends to zero of $ 9 (t) should yield a short range decay. The expression 
in eq. (f3T)|) (given that 7 = 2(f) — 1) tends to a function that scales as £ -1 . This should be 
viewed as a function that scales as S(t) at large tfs, or a short range function. Actually, a 
direct inspection on the right hand side of eq. (f2T)J) recovers this. Since the denominator 
of the right hand side contains the T function, the whole expression vanish as tends to 
zero. Checking more carefully, for = the right hand side of eq. (|2*^|) is proportional to 
exp [— (jJ q t\. If we try now a solution Q q (t) oc exp [— u q t\ we find that it does not work. The 
reason for that and how to obtain the correct asymptotic behavior (eq. E2J)) is detailed in 
ref. Q 

To complete the picture given so far, that is after finding the power law that governs 
the of the structure factor, we need to know the steady state structure factor <ft q and its 
associated "decay rate" uo q that will characterize the short time decay as well. 



B. Steady state properties 

In this part we obtain another equation, that together with eq. (J3*3j) . forms a complete 
set of coupled equations that will yield the small q dependence of 4> q and uj q . In order to 
achieve this we would like to discuss eq. ()27j) in the limit of short-times as well. Here, it is 
more convenient to treat the time-dependent equation directly, so we Fourier transform eq. 
flUj) to yield 

00 

%(t) + {y q -u q ) J dt'e-^'^it-t')- 

—00 

J d d e \A e , q _ e \ 2 J df-^—§ q ^ {t - f) $ £ (t - f) + 



(27T) d 



-00 

00 00 



+ / dHA^A^ IJdt'J dfe-^'-^"^ (f) % (t-t'~ t")+ 



.0 



(2tt) c 

+ J dt' J dfe^'+^'^i (f) % (t-f- t") = — ^ J ^^du. (35) 



-00 —00 



Setting t = and following the same steps described above for long-times (i.e. breaking 
the ^-integration into large and small i regions, and discussing the small q behavior of 
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each) gives the following short-time evaluation of eq. (|35|1 

rqo ,&„ A lA-lKl-<l 



(27T 

2 



(2k) 



E 1( f + J d a 
1 , 

- / \A iA _, 



«*S ({/} 



Ct>a 



COS (7T£ 



(36) 



where as before £i, and -E3 are (renormalization) constants. In addition, we used the 
following notations 



F2 ({/} 



00 00 



Uq-l LOq 
LU£ ' Ld£ 



dx dye—yf (^fx)f(^x + y), 







and 



^3 ({/} 



Uq — l U)£ 
IJJa ' <JJa 



7 (^) / dx - 



(37) 



(38) 



Up to this point we obtained two coupled equations for <p q and u q . Note that the equations 
above depend on the (unknown) functional form of the scaling function /. We will proceed 
now as far as possible without specifying that form, to obtain many general results about 
the exponents. In the actual numerical calculation of the exponents we will resort to an 
approximate form of / to be described later. 

We would like now to solve eqs. (J3*3j) and (}3l)j) in the limit of small c/'s. For convenience 
we rewrite these equations using the explicit form of A(_ m as 



and 



where 



(2v + D l )q 2 -uo q + J< (q) 



(v + E x ) q 2 <\> q + I< (q) <\> q - I< (q) - E 2 



A°° cos (vr0) T (20) ' 



E 3 D° (q) fa)' 2 * 



40-1 



COS KC 



J K (q) 
A< (?) = 



i\ 2 flo 



(2n) d 
3A 2 rio 



d d £ 



?• {q-t)\ [q- (q-£) 



({/} , ^ 



(2k) 



d a 



?• ($-*)] [q-(q-i) 



<>F 2 ({/} , 



and 



(q) 



2A 2 rio 



(2 



K 



d d 



£■ (q-i) 



({/} 



(39) 

(40) 

(41) 
(42) 

(43) 
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Notice that the integrals in eqs. (p£T]) -(pf3 j) are cut by q Q . q is chosen in such a way that 
below it (p q and u) q are expected to be power laws in g, 

cf> q = Aq~ r , (44) 

and 

oo < q = Bq z , (45) 
where z is the dynamic exponent, and T is related to the roughness exponent a by 

a = (r - d)/2. (46) 

As mentioned above, the integrals in eqs. (|41| ) -([43 p are cut by go, and therefore we can 
readily use these power-laws inside the integrals. Using the power laws we can also rewrite 
eqs. (jSnjl - pnj) as 

Aq^ [(2,/ + Di) g 2 - Bq° + J< (q)} = f (g) /^fwl V ^ 

L J A°° COS (7T0J T (20) 

and 

A (i/ + £x) g 2 " r + (g) Ag~ r - 7< (g) - E 2 = ( f^ , (48) 

COS (7T0J 

where we have neglected the E3 term in eq. ()4())1 as it is negligible compared to the left hand 
side in the limit of small g's (since < 1/2). 

It is interesting to notice that these equations are a nontrivial generalization of the Self- 



Consistent Expansion (SCE) developed in refs 



25 



. More specifically, if we take the limit 



of 4> — > 0, and plug in f(u) = e~ u , which is the scaling function of the linear theory when 
= 0, both equations (i.e. eqs. (|4Tj) - (|4*K|) ) reduce to the equations obtained using SCE. 
It is a surprise to find this similarity because the self-consistent expansion was originally 
derived using the Fokker-Planck equation associated with the Langevin-like KPZ equation, 
while the derivation given here deals directly with the Langevin form. Once we realized this 
surprising similarity, it is only natural to follow the asymptotic solution that is used in the 
well-established SCE literature, and is detailed for example in ref. j^j]. 

V. DETAILED ASYMPTOTIC SOLUTION 

As mentioned above, in performing the asymptotic solution of the self-consistent equa- 
tions, we follow previous work. We also focus here, for simplicity, on the case of noise without 
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spatial correlations (i.e. D° (q) = D ). However, eqs. f|^Tj) - f|^H|) are valid for any spatial 
correlations of the noise (i.e. for any D° (g)), so that the more general case is postponed to 
the next section. 

The fist step in the asymptotic solution is to evaluate the integrals I< (q), I 2 (q) and 
J < (q) using the power laws given in eqs. (jSj)-(jlSI) 

. q 2 for c? + 2- r- z>0 

I< (q) , J< (q) oc \ , (49) 

q d +^r-z for d + 2-T - z < 

. const for d + 4 - 2r - 2z (1 - 20) > 

I< (q) oc { V YJ (50) 

?d+ 4- 2 r- 2 y or d + 4 _ 2r _ 2z (1 - 20) < 

We consider now the upper-right quadrant of the (r, z) plane, where a solution may be 
expected. The lines d + 2 — T — z = and d + 4 — 2T — 2z (1 — 20) = divide the quadrant 
into four sectors. We investigate next each sector separately to decide whether a solution of 
the equations ()47j l -(|48 p can exist there or not (in the limit of small q's). 

Sector a is defined by d+2 — T — z > and d + A — 2T — 2z (1 — 20) > 0. In this sector 
equations (J4Tj) and (|4*Sj) reduce to 

and 

A(u + E 1 + Ek) q 2 ~ r -E 5 -E 2 = ^-—^ ; - (52) 

COS (7T0) 

First, the possibility that 2 — T > —2zcj) can be ruled out immediately, because B is 
positive, so that eq. (JSTj) cannot be balanced in leading order (in powers of q). If 2 — Y < 
—2z(p then the right-hand side of eq. ()51|) is negligible compared to the left-hand side, so 
that the leading order equations are identical to those obtained for the white-noise KPZ 
problem and thus the standard KPZ results from refs. 0, 2^ are restored. Therefore, we 
get T = 2 and z = 2. Since T = 2 and z must be positive, the condition 2 — T < —2zd> 

n 

can be met only for < 0. Namely, for the case of noise anticorrelated in time 27J. Such a 
solution holds only for d > 4(1 — 20). 

The other relevant option is 2 — T = —2z<p. This implies that must be positive, T > 2 
and z > 2. There is now an interesting difference between the case z > 2 and z = 2. For 
z > 2, the leading order terms in eqs. (|5T|) and (|52jl lead to two linear homogeneous equations 
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in the quantities A and B This implies that in order to have a physical solution with 
A, B > 0, we must have the determinant of the coefficient matrix vanish, namely 

T (20) A°° {2v + D x + D 2 ) = (u + E 1 + E A ) , (53) 

Since the quantities Di, D 2 , Ei, E 2 depend on the behavior of (pi and u>i for £ > q , on 
the total upper cutoff etc., it is difficult to envisage that eq. ()53j) can be fulfilled under 
accidentally, for non-generic values of the parameters of the system. The case with z = 2 
is different. The two equations for the coefficients A and B have now an additional term - 
AB on the right hand side of the first equation. This enables now a generic solution for the 
coefficients. In that case T = 2 + 40. Considering the defining conditions for the sector we 
find that within the sector such a solution is possible only for d > 4. 

Sector (3 is defined by d + 2 - T - z > and d + 4 - 2r - 2z (I - 20) < 0. In this sector 
equation (JH2j) is replaced by 

7~) f) — 2(f> 

A(u + E 1 + E<) q 2 ~ T - £ 6 g d+4 - 2r -* - E 2 = ^hr^ 2 , (54) 

COS (7T0J 

while equation (|51|) remains intact. 

The analysis for the possibility 2 — T = —2z<j) in sector (3 is similar to the above analysis 
for sector a. The only difference is that due to the different defining conditions of the sector, 
such a solution with z = 2 and T = 2 + 40 holds within the sector for 2 + 40 < d < 4. 

Combining the results for sectors a and (3, we see that for < the noise term is 
irrelevant, and the critical exponents that describe the EW problem with uncorrelated noise 
(namely T = z = 2) are restored. This option is possible if d > 2. In addition, for > 
we get the new solution z = 2 and T = 2 + 40 that is just the solution obtained for the 
EW equation with temporally correlated noise (see eq. ()13|) above). Following the above 
discussion it is realized that this solution is possible only for d > 2 + 40. Therefore, the 
lower critical dimension in this problem is d c = 2 + 40 (provided > 0, otherwise d c = 2 as 
mentioned above). 

Sector 7 is defined by d + 2 — T — z < and d + 4 — 2T — 2z (I — 20) > 0. In this sector 
equation is replaced by 

7") f) — 2(f> 

A(v + Ek) q 2 - T + AE 7 q d+4 - 2r - z - E 5 - E 2 = ^-^q-^. (55) 

COS (7T0) 

First, The two defining conditions of this sector imply that the first term on the left hand 
side is negligible compared to the second term, and the second term is negligible compared 
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to the term on the right hand side of the equation. Therefore, looking at the simplified 
equation, we must conclude that = and — E 5 — E 2 = . However, this is impossible 

because the left-hand side is negative definite. 

Sector 5 is defined by<i + 2- r- 2<0 and d + 4 - 2V - 2z (1 - 20) < 0. In this sector 
eqs. (JUj) and pKj) take the form 



Aq- r 



> + D x ) q 2 - Bq z + 



8A 2 A 



d+4-r- 



z G({f},T,z) 



D B- 2 ^ 



cos (vr0) r (20) 



-2<t>z 



and 



A (v + q 2 - T - J^g^-'F ({/}, T, z) - E 2 " 



DnB- 2 * 



(2n) d B 
where G ({/}, T, 2) is given by 



G({f},T,z)=J d d t 
and F({/},r,z) given by 
F({f},V,z 



COS 7TC 



:9 



-2<t>z 



(56) 
(57) 









e- (e-t) 




e-t 


z 



e - 1 



-4-r 



(58) 



+ 



I / cf*t 
eft 



F- (e-t) [e- (e-t) 



e-t 



-t-r 



e-t 



>3 ({/} 



Mm 



e-t 



(59) 



e is a unit vector in an arbitrary direction, and the t-integration is over all (^-dimensional 
space. 

From the defining conditions of this sector, it is possible to neglect the g 2 -term in the 
brackets on the left-hand side of eq. (|56j) compared to the third g d+4 ~ r ~^-term (since d + 
4 — T — z < 2). In addition, It is also possible to neglect the g 2_r -term on the left-hand side 
of eq. (|H7jl compared to the second g d+4_2r_2 -term (since d + 4 — 2T — z < 2 — T). It is easy 
to see that for < all the usual KPZ results are trivially retained, as the right hand side 
is irrelevant then. Therefore, we focus on the case > and thus the constant term on eq. 
(15 7|) can be neglected compared to the noise term on the right hand side. 

Now there are two options: First, when r > z (1 + 20) then the right-hand side of both 
equations (jHEJ)-(jSZ|) is negligible compared to the left-hand side. In that case the critical 
exponents are determined by a combination of the scaling relation d + A — T — 2z = 0, and 
the equation 

F[{/},r,z(r)] = o, (60) 
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(where F is given by eq. (|59~j) above). We denote the solutions of the transcendental equation 
by T ( f ) (d) (since the exponent T is dependent on the spatial dimension d and on 0). For 
example, in one dimension, and for = it can be shown analytically that Tq (1) = 2 and 
in two dimensions a numerical solution of the equation (again for = 0) yields Tq (2) = 2.59 
(see refs. j^ll^). A discussion for general 0's will be given below. Still, we must remember 
that a solution here is obtained by requiring r > z (1 + 20). This yields a necessary condition 
for the existence of such a solution, < ,; 



2(d+4-rv(d)) ' 

The second option in this sector is Y = z (1 + 20) (the possibility T < z(l + 20) is 
irrelevant because one cannot balance the equations and still be consistent with the defining 
conditions of sector 5 in that situation). Then in order to balance eqs. l)56J) - (|57jl we must 
also have d + A — 2Y — z = —2z(f). This leads to the new solution z = and Y = 

(d + 4) However, this solution is valid only if the equation (JoUj) does not yield exponents 
(namely Y^ and z^) that make the g d+4_2r_z -term in eq. ([57]) dominant. Not surprisingly, 
this requirement translates into the condition > ^{dt^ ? ^i)) ~ i m Ply m § either a smooth 
transition between the two types of solutions, or a complete domination of the first option 
(i.e. r^(d)). Actually, the existence of this new solution also requires F [{/}, Y, z (Y)] < 0. 
This requirement turns out to be the same as Y^d) < Y new = (d + A) |^f|, so that this 
extra requirement is fulfilled automatically since Y = z (1 + 20). 

To summarize the results of sector 5, we found two possible strong-coupling solutions. 
The first solution is obtained from the equation (}60|) . and its scaling exponents are denoted 
by Y<f,(d) and z^d) (this solution reduces to the standard KPZ results when = 0). The 
second solution, is given by the explicit expressions Y = (d + 4) and z = j^;. Then, 
in a given dimension d and for a given 0, the actual strong coupling exponents of the KPZ 
problem with temporally correlated noise are just Y = max{Y ^(d), (d + 4) |-|f^} and its 
corresponding z. Thus, the transition between the two solutions as a function of (if such a 
transition exists) is continuous. However, it should be emphasized that for a specific one 
of these solutions dominates so there is no phase transition between them. 

Based on the r e sults „ t th e lln e ar th e 0ry ^ and the nn m e rlC a, slm n,a tl on Q we expect 
the exponent T to be a nondecreasing function of (that means that the inclusion of temporal 
correlations does not make the surface smoother). This implies that Y^d) > Y (d), and since 
we know Y Q (d) from the white- noise KPZ problem it is easy to determine a lower bound on 
0, denoted by c (rf), such that the new solution Y = (d + 4) is not possible below it. 
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The specific value of (j) c (d) is 4> c (d) = 2 ^^ 4 d |~^^ . The fact that T^(d) is non-decreasing as 
a function of implies two possible options. Either the expression T^d) gives the strong 
coupling solution for the whole range of < < 1/2, or it is the solution for small and 
crosses over to V = (d + 4) Such a crossover can occur only above <fi c (d). 

We turn now to the actual evaluation of the exponents. To do that we need an ansatz 
for the scaling function / on which the form of the equations for T and z depend. Since the 
equations were constructed in such a way that second order corrections to the quantities 
<p q and u q vanish, we use as in refs. j^, [^J , the zero order form of the scaling function in 
evaluating these corrections. The scaling function in zero order is the function obtained for 
the corresponding linear theory, given by eq. 1)170 . We therefore simplify eq. (J5Uj) using this 
ansatz. First, one quantity can be evaluated exactly 

Fi ({f E w} , a) = / e~ ax f E w (x) dx = - - . (61) 
jo cr — 1 

In addition, the functional F 2 can be simplified so that it involves only one dimensional 
integration (instead of double integration) 

00 00 

F 2 ({jW} , a, b) = J dx J dye~ x ~ y f E w (ax) f E w (bx + y) 

du^-—- — - ^r 2 Mi-*H (62) 





cos (nd 



| u |-2* x a l+20( 


1 - ibuY 24, - (1 


— ibu) 


v? + 11 — iu 


a? — (1 — ibu) 2 





(note that the integral is real, as it should be, even though the integrand is complex). 

In order to illustrate the outcome of this analysis we specialize to one dimension. First, in 
one-dimension only strong coupling solutions are possible as the critical dimension is 2 + 40. 
Second, as mentioned above, in one dimension eq. ()6T])l can be solved analytically for = 
and it gives r = 2. This corresponds to a roughness exponent of cto = 1/2 (using eq. flUH) ) 
and to a dynamic exponent of Zq = 3/2 (using the scaling relation z (T) = (d + 4 — T)/2). 
for higher values of <fi one has to solve equation (jBHj) numerically using the ansatz of the 
linear theory. These results are summarized in Fig. El as the solid line. The figure also 
presents the possible second solution Y = 5(1 + 20)/(3 + 20) (that corresponds to a = 
(1 + 40)/ (3 + 20)) and z = 5/(3 + 20), and a continuation of r as a dashed line. However, 
since the this solution is smaller than r^, it is practically irrelevant, since dominates the 
whole range. 
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FIG. 2: (a) The roughness exponent a^, (b) the growth exponent (5^ and (c) the dynamic exponent 
as a function of the exponent <j> for decay of temporal correlations in d = 1. Note that the 
dynamic exponent z was inferred for the numerical results of ref. ^3] from a and [3 using the scaling 
relation z = a//3. Second, notice that the DRG result is possible only up to (j> = 0.46. Third, the 
dashed line shows our second possible solution (using SCE) that turns out to be irrelevant here, 
since it is smaller than the SCE for all </>'s. 

It is particularly interesting to compare this one- dimensional result to the DRG result 
that was presented in the introduction [lOj, and plotted for convenience in Fig. Generally 
speaking, the two methods disagree on the values of the critical exponents significantly 
over most of the 0-range. Three substantial differences can be observed between the two 
methods. First, using the self-consistent approach we found no "threshold behavior". That 
is, we found a continuous variation of the scaling exponents a and z as a function of (f> over 
the whole range of possible 0's, rather than no variation of these exponents up to a critical 
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value of 4> c and a quasi-linear behavior from that point on. 

Second, we found a solution for the exponents for every <p, while the DRG approach 
found no solution above <fi = 0.46 (claiming that no stable surface can grow under the 
condition 0.46 < <p < 0.5). Interestingly, the threshold <p c (the crossover point) that was 
predicted using DRG in d = 1 (namely 0.167) is the same as the lower bound we found above 
(0 C = 1/6 for d = 1). Our exact statement was that for <fi > 1/6 the second strong-coupling 
KPZ solution becomes possible in principle (but not in practice). Therefore, the DRG result 
might reflect this exact statement. 

Third, we found that z is a decreasing function of 0, while the DRG approach predicts 
an increasing value of z. The reason for this difference is not clear, but it might stem from 
the definition of the typical "decay rate", that was actually defined using the scaling form 
(124)1 . rather than a more "traditional" definition such as eq. (119)1 . The reason for using this 
definition is that the integral over the scaling function f(u) does not converge, because of its 
power-law tail. Actually, in the case of the linear theory, where everything can be calculated 
exactly, the only possible definition is the one we used. Now, since the introduction of 
temporally correlated noise certainly slows down the relaxations in the system, this might 
have caused an artifact of increasing z, because larger z's are interpreted as longer relaxation 
times. However, in our approach, we do see this slowing down clearly, but it does not come 
from a larger dynamic exponent z in an exponential decaying scaling function, but rather 
from a very slowly decaying scaling function, which does not decay exponentially. Thus, 
this difference might reflect a better understanding of the time-dependent dynamics in such 
driven systems. 

The discrepancy between DRG and our result is obvious. This should not come as a 
surprise as it is not new that SCE and DRG yield different results, apart from some special 
cases, like the one-dimensional white noise (in space and in time) case. In fact it should 
not be worrying in the present case because in many cases studied in the past the results of 

ations 



SCE are either dimilar of those of DRG (e.g., the one dimensional case with no corre 
QQ), or superior to them (as in the case of white noise KPZ in higher dimensions 



I , or 

3)- 



in the case of the nonlocal KPZ equation that is exactly soluble in one dimension 
On that basis we may expect this to be true also here. This expectation is supported, in 
fact, by the non physical feature of increasing z predicted by DRG (see Fig. |2J). 

We must admit that we were surprised by the fact that that our results do not show 



22 



a threshold behavior (like DRG). Namely, the fact that we do not find the characteristic 
exponents for a range of small 0's to be identical to those found at = 0. In fact, it 
is interesting that such a solution to the SCE equations exists but it is only potentially 
possible, the actual solution that determines the exponents is the second solution (in one 
dimension). The characteristic exponents obtained by the scaling approach [18] and by the 
numerical simulations do not exhibit a threshold behavior as well (see Fig. |2J). A closer 
inspections of the values of the results of the scaling approach, reveals that these results are 
actually a simple interpolation (linear for f3 and almost linear for a and z) of the values of 
the exponents between the results we predict for = and = 1/2. 

The results of the simulations are even more interesting. The simulations jl2j reports the 
values of a and (3 as a function of (these values are recovered in Fig. EJ). The first four 
small-0 points for a and the first three small-0 points for (3 agree with our curves. Thus, 
the simulations and SCE, which are completely independent, agree exactly in the region 
of where threshold behavior could have been expected. This suggest that the small 
behavior predicted by both methods is correct and indeed a threshold behavior should not 
be expected here. For higher 0's the results of simulations obviously deviate from our curves. 



Can the results of simulations be trusted for large 0's? In their paper Lam and Sander [17 1 
report difficulties observed for larger 0's, but they also report measures taken to ensure 
the correctness of their final result. To check whether indeed those larger results could 
be trusted, we took the freedom of drawing the dynamical exponent z = a//3 as inferred 
from the simulations that give only a and (3 (the inferred values of z are presented in Fig. 
Efc) above). The erratic oscillations of the dynamic exponent obtained from the simulation, 
strongly suggest that the measures taken by the authors to eliminate the observed larger 
problems were probably not enough. 



VI. GENERALIZATION TO SPATIO-TEMPORALLY CORRELATED NOISE 

Up to this point we discussed the relatively simple case of noise without any spatial 
correlations (i.e. D° (q) = D ). However, including spatial correlations bears no principal 
difficulty to the analysis presented above. For example one can just replace Do with D° (q) = 
-Do<?~ 2p from eq. (|51jl and on, and thus can easily obtain the scaling exponents for that case. 
For the sake of presenting a complete picture we briefly summarize the results obtained for 
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this case. First, there is the weak-coupling solution, which is again just the corresponding 
EW result for such a noise, given by 

z = 2 and V = 2 + 2p + 40. (63) 

The weak-coupling solution is possible for d > 2 + 2p + 40, so that here the lower critical 
dimension is d c = 2 + 2p + 40. 

Second, there is the strong coupling solution, given by 

**(<0 d V= [ v +W 2 P + (i + 20) H (d) < r (d) 

^sr an 1 (d+4 - 3 s i+20) 2 P+ (i+20)^(d)>r,(d) ' 

where as before, (d) is the solution of the equation (|6T)|) and (d) — (d + 4 — (d))/2. 

Furthermore, the method presented above is not restricted to noise terms that have sep- 
arable correlators, i.e. D (q, u) = D° (q) uj~ 2< ^, and can just as well deal with non-separable 
correlators (that is any functional form of D (q,u)). In that case, the only difference would 
be to replace the right-hand side of eq. (|2*7j) by the expression ^ 2 +'^2 with the required 
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D (q, lu) inside. 

In order to demonstrate this option, we discuss an interesting application of this approach 
to the KPZ equation with a very special kind of spatio-temporally correlated noise (this result 
was mentioned at the end of the introduction). This problem was previously solved by Li et 
al. Q| in the context of Vortex lines in the three-dimensional XY model with random phase 
shifts, and it boils down to solving the KPZ equation in two dimensions (more specifically 
the 2 + 1 case) with a noise term that has the following spatio-temporal correlations 

( V (f, t) V (?, 0> = , 2 , (65) 

where r is a two dimensional vector (i.e. f = (x,y) G R 2 ). In order to apply the method 
presented above for finding the critical exponents of this model we have to Fourier transform 
the noise correlator first 

x2 5 2 (q + <f)6(u + u') 
q 2 + uj 2 

(where qis a two-dimensional vector) so that D (q,u) = 2+^,2 (with Do = a (Jm) 2 ). Now, 
for small q's u q = Bq z , and we can evaluate the Fourier integral explicitly (assuming z > 1) 

iut 1 7T e~ u ^ 1 ire^ , , 

duJ ~ 7^0 -ttt = 00 • ( 67 ) 



(7](q,uj)r](q ,u)) = a (Jm) — — , (66) 



(u 2 + u 2 ) {uj 2 + q 2 ) ' u$2q*/u* <? 2+z 25 
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As mentioned above, this expression should replace the right-hand side of eqs. (|29|) and 
([35)1 . One can easily be convinced that this term is subdominant in the long-time limit, and 
therefore drops out of eq. (jSHJ)- However, in the short-time limit this term is not negligible, 
and therefore modifies eq. (jSfij) accordingly. 

Following all the required steps from that point on leads to the results T = z + 2 and 
d + 4 — 2z — T = 0, so that for the case we were interested in (d = 2) we get 

z — — and T = — , (68) 



identical to the analytic result presented in ref. 19] (a numerical analysis presented there 
also verifies this result). 



VII. SUMMARY AND CONCLUSIONS 



In this paper we developed a time-dependent self-consistent approach to deal with the 
KPZ equation driven by a temporally correlated noise. This achievement was made possible 
thanks to the observation that there is a time scale separation between short-time and 
long-time behavior of the system. More specifically, it was realized that when temporally 
correlated noise is present in the system, then slow relaxations of various time-dependent 
quantities should control the long time behavior (in this case algebraic decay of the time- 
dependent correlation function $ g (t)). In addition, it was seen that the short time behavior 
is influenced by the long time behavior and vise versa. 

To summarize the results briefly, we found that the KPZ equation with temporally cor- 
related noise, just like the problem with white noise, has both a strong-coupling and a 
weak-coupling solution. The weak coupling solution is described by the scaling exponents of 
the corresponding linear theory (EW equation), and is made possible for dimensions higher 
than the lower critical dimension d c = 2 + 40 (the specific values are given in eq. (JBlIj) . The 
strong-coupling solution, which is relevant also for low dimensions, is described by critical 
exponents that are a result of a competition between two possible solutions. First, there 
is an extension of the classical white- noise KPZ solution denoted by T^d) that is derived 
from an integral equation ([60)1 . The other possible solution is a new strong-coupling solution 
and is given in eq. ([64)1 . The actual exponent T that describes the surface is the maximum 
between the two options. For small 0's T^(d) is the solution but it may (or may not, as in 
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the one dimensional case) cross over in a continuous manner to the second solution. For a 
detailed discussion and comparison to other methods in one dimension see section V. 

The comparison suggest that for small values of our results are correct, being supported 
by the totally independent numerical simulations. For larger values of 0, the erratic behavior 
of z obtained from the simulations suggests a problem in the numerical evaluation of the 
scaling exponents. Therefore, an independent study, preferably a heavy numerical study, 
aimed at larger 0's could prove useful in a critical determination of the scaling exponents. 
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